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The quantum nonlinear dimer consisting of an electron shuttling between the two sites and in 
weak interaction with vibrations, is studied numerically under the application of a DC electric 
field. A field-induced resonance phenomenon between the vibrations and the electronic oscillations 
is found to influence the electronic transport greatly. For initially derealization of the electron, 
the resonance has the effect of a dramatic increase in the transport. Nonlinear frequency mixing 
is identified as the main mechanism that influences transport. A characterization of the frequency 
spectrum is also presented. 
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INTRODUCTION 

The response of an electron-phonon system to an ex- 
ternally applied electric field forms an important area of 
research in condensed matter physics [1—5]. For exam- 
ple, it is known that coupling between the electron and 
phonons leads to derealization of the initially localized 
Wannier-Stark states under the application of a constant 
electric field (see for e.g. [6, 7]), and a nontrivial response 
to a sinusoidally varying field applied to an electron in 
a lattice results in dynamic localization of the electron 
[8-11]. By tuning the Wannier-Stark splitting to overlap 
with the lattice vibration frequency, new optical phonon 
mediated paths have been observed experimentally in su- 
perlattices [12]. Advent of novel terahertz spectroscopy 
techniques has made it possible to observe nonlinear re- 
sponse of systems to applied electric fields in real time 
along with the small energy scales. Experimental obser- 
vation of the evolution of long-range correlations leading 
to the birth of the quasiparticles [13] , resolution of the dy- 
namics of the nonlinear transport of polarons in the pres- 
ence of strong quasi-stationary fields [14], ultrafast res- 
olution of the elementary excitations in superconductors 
[15], observation of the strong mode-coupling dynamics 
in the cavity quantum electrodynamics [16] for example, 
are some of the recent hallmark demonstrations of the 
modern experimental sophistication to resolve real-time 
response of the modes of the condensed matter systems 
to the external perturbations (e.g. electric fields). Moti- 
vated by these recent experimental progress, this paper 
numerically investigates electronic transport in a simple 
dimer system, consisting of an electron oscillating be- 
tween two sites in the presence of the applied DC electric 
field. The electron, described quantum-mechanically, is 
allowed to couple weakly with the vibrational degree of 
freedom, represented by a classical simple harmonic os- 
cillator. Henceforth in this paper we will refer to the 
oscillations of this simple oscillator as a phonon field or 
lattice vibrations. 



Most generally, the system is described by a Hamilto- 
nian that separates the electron phonon system of inter- 
est into three parts [17, 18], 

Htot = H e i + H p h + Hi n t (1) 

where the first term describes a tight-binding electron in 
a 1-dimensional lattice: 

H el = \m)(m + 1| + \m){m - 1| + E S^m\m){m\- 

m m 

(2) 

with the lattice site integer m running from — oo to +oo, 
V is the nearest neighbor transfer rate, and E = ea£ is 
proportional to the magnitude of the degeneracy lifting 
electric field £ with proportionality constant as a prod- 
uct of electronic charge e and a lattice constant a. The 
phonon term in Eq. (1) is represented by a collection of 
identical Einstein oscillators of mass M and frequency 
w , 




m m 



where, x m and p rn are respectively, the position and mo- 
mentum of the mth oscillator. The interaction (Eq. (1)) 
between the electron and the lattice is taken to be linear 
in the displacement of the oscillator coordinate and to 
modulate only the site energy of the electron and not the 
intersite transfer matrix elements: 

Hint = ce*^2 \ m )( m \ x m- (4) 

m 

The parameter a describes the strength of the interac- 
tion. 

In this paper we focus on the problem in which the 
lattice is restricted to only 2 sites. The 2-site system, 
also known as the spin-boson system or the 2-site Hol- 
stein polaron [19, 20], has been explored extensively due 
to its formal simplicity and capability to reproduce cer- 
tain features observed in extended systems rather accu- 
rately [21, 22]. Despite the simplicity of the system, the 
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apparent complexity in its behavior is evident from the 
fact that the complete quantum mechanical problem can 
not be solved exactly analytically. Perturbation methods 
and approximation schemes are thus commonly used to 
obtain insight into the system. In the present work we 
concentrate on the semiclassical approximation scheme 
[23-27], the range of validity for which was investigated 
elsewhere [28, 29]. 

In the semiclassical approximation, where only the 
electron is treated quantum-mechanically, the equations 
of motion are derived by assuming for the electronic 
states \ip e ){t) = J2 m c m(t)\m), where c m (i) is the prob- 
ability amplitude to find the electron in the localized 
state |m). The time evolution of the probability am- 
plitudes c m (t) are obtained from the Schrodinger equa- 
tion (ih d\tjj e )/dt = {H e i + H in t)\tp e })- For the equation 
of motion for the lattice one uses the classical Hamilto- 
nian (H) = (?/' e |-fftot|V'e), and the Hamilton's equations 
x n = d(H)/dp n andp„ = —d(H)/dx n . In such case, one 
obtains, 

ihc m = V (c m+ i + c m _i) + ax m c m 

M 

The equations of motion governing the motion of the elec- 
tron and lattice for a 2-site system are immediately ob- 
tained from Eqs. (5) and are given by, 

ihci = Vc2 + ax\C\ — Eci (6a) 

ihc 2 = Vc\ + ax 2 c 2 + Ec 2 (6b) 

Xi +ul Xl = -|||ci| 2 (6c) 

x 2 +uj 2 x 2 = -§ |c 2 | 2 (6d) 



Various variants of Eqs. (7) have appeared extensively 
in the literature. For example, under the adiabatic ap- 
proximation (y = p), Eqs. (7) can be written in a closed 
form in p as, 

p+(l + A 2 + XC0 ) P = -A C0 -^ 2 -^ P 3 (8) 

2 

with Co = f*o — ^y 1 — Apo as a constant which depends 
only on the initial conditions and the system parame- 
ters. Eq. (8) is similar to the equation for a degenerate 
trimer the exact solution to which in the form of Weier- 
strass Elliptic functions has already been obtained [31]. 
Setting A = in Eq. (8) reduces it to the well known 
Discrete Nonlinear Schrodinger equation (DNLSE), for 
which the exact solution is in the form of Elliptic func- 
tions [21]. Amongst other important studies the response 
to different initial conditions [32, 33], calculation of line 
shape in neutron -scattering [34] etc. have been studied. 
Attempts have also been made to address the issue of 
validity of the adiabatic approximation for a degenerate 
dimer. Eliminating the approximation renders the prob- 
lem extremely complicated, making it difficult to obtain 
analytical solutions. In certain instances, some analytical 
progress has been made either by introducing dissipation 
[22], or obtaining exact solutions for the non dissipative 
case, valid, for a limited set of initial conditions and pa- 
rameter values [30]. Our goal in this paper is two-fold: 
firstly, we do allow the lattice to vibrate on its own time 
scale and secondly, we focus on the response of the elec- 
tron in the non adiabatic dimer to the application of elec- 
tric field concentrating specifically on its transport. 



Pauli spin matrices a x , a y , & z here p, q, r respec- 
tively, can be obtained through p = en — C12, q = 
i(ci2 - c 2 i), r — C12 + c 2 i with c y = Cjcj. Fur- 
thermore by defining a dimensionless internal coordinate 

V = ( x 2 — xi) Eqs. (6) can be recast in terms of real 
quantities only to give [20, 22, 30], 

p = q (7a) 

q = —p — xu r — Ar (7b) 

r = xyq + Ag (7c) 

y + u 2 y = uj 2 p- (7d) 

where p implies a derivative with respect to a dimension- 
less time r = (2V/H) t. The various dimensionless param- 
eters appearing in Eqs. (7) are \ — « 2 / (2Muj 2 V) , A = 
E/V and cj = Hujq/(2V). The main observable — 1 < p < 
1 describes the difference in the probability of occupancy 
of the two sites of the dimer, while, y is the scaled in- 
ternal coordinate of the lattice displacement. The limit 
of weak coupling of electron to the lattice corresponds 
to x ^ 1 ancl we w iH restrict our analysis only to this 
regime. 



RESULTS 

To quantify electronic transport in the dimer model 
of interest, it is important to identify time- independent 
(steady-state) parameters describing the dynamics of the 
system. The importance of time averaged probability dif- 
ference (p) = lini^oo 1/T J Q T g?t p(r) has been realized 
previously [31]. While (p) is a useful quantity in describ- 
ing characteristics of the dynamics, it does not generally 
yield information about the transport without additional 
knowledge of the magnitude of the electronic oscillations. 
In this spirit, we introduce a quantity given by 

r (i-fw/<p», (p)<o, 

Id = < oo, (p) = 0, (9) 

I 1 Pmin /(p), (p>>o, 

where p ma x(min) is the maximum (minimum) value of 
the probability difference p. Note that < 7^ < 1 im- 
plies localized dynamics while 7^ — > 00 together with 
Pmax(min) 7^ implies delocalized dynamics. We there- 
fore refer to as the degree of transport, which is also 
characteristic of the electrons mobility in this context. 
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To demonstrate the usefulness of the 7,2 measure, let 
us first turn off the lattice coupling by setting \ = in 
Eq. (8). Application of the electric field beyond a certain 
value would lead to the localization of the electron in the 
Wannier Stark states, and hence, to hindered electronic 
transport. This can be understood simply by considering 
a constant energy picture obtained from Eq. (8) 

■ 2 

V — + U(p) = E constant (10) 

2 

where potential energy U (p) — ^fp 2 + Acop, with cu e = 
\/\ + A 2 being the bare-electron frequency. For asym- 
metric initial condition (po — 1), (p) and p m in can be 
trivially calculated from Eq. (10) to be 1 ^ A2 and ^2+1 
respectively, which when substituted in Eq. (9) gives 
7d = i. Thus, in the Stark case lies in the range 
[0,oo) depending upon the field strength: for A — » 0, 
7d — > oo showing complete transfer of the particle on 
both lattice sites with the probability difference oscil- 
lating between the two possible extrema ±1. Applica- 
tion of field results in a nonlinear reduction of the de- 
gree of transport with increasing field. A = 1 results 
in the localization of the electron on one of the lattice 
sites. Further increase in field value reduces the am- 
plitude of oscillation and is effectively captured by jd- 
In the other extreme if we start with completely de- 
localized initial condition p — 0, r = 1 which gives 
(p) = ~j^2 and p m ax = 0, for which 7<j = 1 indepen- 
dent of the value of A. While 7^ correctly predicts inhib- 
ited transport (particle now oscillates only on one site), 
it does not capture the variation in the amplitude of os- 
cillation as it is insensitive to the magnitude of the ap- 
plied field. Similarly, in the degenerate adiabatic non 
linear dimer initially localized conditions results in the 
trapping of the electron on one lattice site as coupling 
parameter x equals 2 [21]. Further increase in field val- 
ues lead to reduction in amplitude of the oscillation and 

is also captured by 7,2 which equals x ^ x 4 for x > 2. 

The other extreme with tq = 1 is a stationary solution of 
the problem leading to no oscillation in p. Breaking the 
degeneracy for this initial condition results in the con- 
finement of the particle on single lattice site. It should 
be noted that if the lattice coordinate is not assigned 
any degree of freedom then completely delocalized initial 
condition leads to either confinement on one lattice site 
(A > 0, x = 0), or stationary solution (A = 0, x > 0) 
involving no oscillations in p. Even for A > 0, x > 0: this 
initial condition simply leads to reduction in amplitude of 
oscillation apart from confinement. This should also be 
evident from Eq. (8) which in the potential picture rep- 
resents the motion of a particle in a fictitious potential 

given by U(p) = Ac p+ ^^^p 2 + ^p 3 + ^p 4 . In- 
crease in x or A has the effect of increasing the potential 
barrier for the electron. 



We now proceed to analyze the effect that applied elec- 
tric field has on a quantum dimer with the lattice allowed 
to vibrate at its own characteristic frequency. For this, 
we consider the case of delocalized electronic initial con- 
ditions (ro = 1) which corresponds to a trapped state in 
the adiabatic case and a stationary state for the degener- 
ate non-adiabatic case. We keep the lattice to be initially 
unexcited (y = 0, y = 0) and stay in a regime in which 
the coupling between the lattice and the electron is weak 

(X« !)■ 
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FIG. 1: A surface plot of the mobility of the quantum particle 
7,j(A,tj) for x — 0.005 and initially delocalized conditions 
ro = IjPo = Qo = Vo = Vo = 0. A sharp increase in mobility 
(white region) is observed along y/l + A 2 ~ uj. 

We investigate the behavior of the system numerically 
by analyzing the dependence of degree of transport 7,2 
with varying electric field A and lattice frequency uj. In 
Fig. 1 we show a contour plot of 7,2 (A, a;) via simula- 
tions of Eqs. (7) for the above mentioned initial condi- 
tions. First, given field value A (y-axis) determines the 
bare-electron frequency u 2 = 1 + A 2 . By scanning the 
lattice frequency uj (x-axis), the behavior in adiabatic 
(uj e <C uj) and uncoupled-lattice (uj e uj) regimes can 
be observed. Indeed, both regions show inhibited trans- 
port as jd = l(Fig. 1, black color), consistent with the 
corresponding analytic results discussed above. Due to 
the presence of the interaction term x, the behavior of the 
system is however non-trivial when the two frequencies 
are in close proximity (Fig. 1, white color) resulting in 
strong enhancement of the electronic transport. In this 
region both oscillators are resonantly coupled, exhibiting 
coupled mode behavior. In other words, the bare-electron 
frequency gets re-normalized upon being dressed by the 
interaction with the lattice. 

This point is further made clear in Fig. 2 where the 
time evolution of both p (solid line) and y (dashed line) 
is plotted for cases (i) "away" (Fig. 2 (c)) and (ii) "at" 
(Fig. 2 (d)) the condition of resonance, as represented 
by the dashed vertical lines in Fig. 2 (a). While on reso- 
nance, energy of electron motion gets transferred into the 
lattice mode, as shown by substantial increase in lattice 
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FIG. 2: Resonant enhancement in jd (a) and normalized lat- 
tice amplitude A\ (b) with varying electric field A for u — \/2, 
X — 0.005 and same initial conditions as in Fig. 1. Time evo- 
lution of p (solid line) and oscillator coordinate y (dashed 
line, normalized to resonant value in (b)) plotted for: (c) 
A = 0.5 away from resonance (vertical line labeled (i) in (a)); 
(d) A = 1 (label (ii) in (a)) for one oscillation period of the 
difference frequency ojdfg- 



oscillator amplitude Ai (Fig. 2 (b)). Lattice motion in 
turn leads to qualitative change of the electronic motion, 
resulting in un-trapped oscillations of the electron be- 
tween two extrema of probability difference p, and hence 
enhanced transport. 

To understand the nature of the coupling we consider 
the behavior of the system in the frequency domain. 
To this end, we obtain frequency dependent probability 
difference V(uj e ) by taking the real part of the Fourier 
transform V(uj e ) = 1ZT{p{t)\ of the time dependent so- 
lution p(r). In Fig. 3 (right) we show the contour plot of 
V(uj e , A) for a fixed lattice frequency. As has been dis- 
cussed above and is evident from 7d(A) plot (Fig. 3(left)) 
transport is enhanced only in the region of coupling i.e. 
when the two frequencies uj e and uj are in close prox- 
imity. What is evident from the contour plot, however, 
is also the presence of new frequency components which 
result from the nonlinear mixing of two fundamental fre- 
quencies {uj e and u>). New low frequency components 
occur when the two frequencies are at resonance, or in 
other words, the coupled modes are near the anti cross- 
ing point. It will be shown in the next section that the 
low frequency contribution is due to second order non- 
linearity, resulting in the difference frequency generation 
(DFG): wdfg = w e ~ which approaches but never 
touches zero due to the anti crossing nature of the cou- 
pled modes. It is this difference frequency term that is 
responsible for the enhancement of the transport on res- 
onance, as this is the frequency with which electron tun- 
nels between the dimer sites, assisted by the phonon field. 
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FIG. 3: Spectral response of the probability difference as a 
function of the electric field (right) . The value of the coupling 
X is 0.005 while the vibrational frequency of the lattice is fixed 
at v2. On the left is shown the degree of transport (jd) with 
respect to the field. The dashed line at 7^ = 1 shows the 
absence of transport in the DC Stark case. Other parameters 
are the same as in Fig. 1. 



Characterization of Frequencies 

In the previous section we have seen that as the res- 
onance condition is met, it results in the generation of 
higher order frequency components, in particular, the dif- 
ference frequency component wdfg (see Fig. 3). The new 
frequency components will in general carry signatures of 
the nature of coupling between the lattice and the elec- 
tron present as a result of the linear interaction in the 
Hamiltonian (Eq. (2)) and the semiclassical approxima- 
tion resulting in Eqs. (7). To characterize the different 
frequency components present we first notice that start- 
ing from Eqs. (7) two integrals of motion can be derived 
[30], 

p 2 +q 2 +r 2 = 1, and (11a) 
v 2 v 2 r A 

±^ + y py+ p = X- (lib) 

2uj 2 2 XX 

where Eq. (11a) simply means that the total probability 
is always conserved, whereas Eq. (lib) is the statement 
of conservation of energy in the system. The presence 
of these two integrals of motion in Eqs. (11), reduces 
the number of independent degrees of freedom to three, 
which still does not permit exact solutions. Solving for r 
in Eq. (lib) and substituting in Eqs. (7), we obtain the 
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following second order equation in p 



p= -p 



l + A 2 {l + ey} 2 -A(l + ey) X I 



(1 + ey) (y 2 + £ 



(12) 



with e as the ratio of interaction strength x an d the field 
strength A. For situations in which the applied field is 
much stronger than the interaction strength we set e = 
in Eq. (12) to obtain, 



P- 



1 



A 2 )p 



Ax/ • \f ( 



V 



(13) 



which for our choice of initial conditions 
{Vo = Vo = 9o = 0) and by substituting 1 = r /x-Ap /x 
reduces to 



1>- 



l+A 2 ) P 



-A(r -Ap )+^(,r 



V 



(14) 



The set of equations given by Eq. (14) and Eq. (7d) are 
able to reproduce the dynamics of the original equations 
(Eqs. (7)) rather accurately for small values of X - Hence, 
in order to obtain an estimate of the different frequency 
components generated, instead of solving Eqs. (7), we 
take Eq. (14) and Eq. (7d) as the point of departure. It 
is clear that the two fundamental frequencies in equations 
(14) and (7d) are given by ui e — vT-j- A 2 for the electron 
and u> for the lattice parameter y. Using a = xA/2 as the 
perturbation parameter we look for asymptotic solutions 
in the form of, 



P 



= E« fe 



(15a) 
(15b) 



where k is the index of the power series expansion and 
takes values 0, 1, 2. .etc. Substituting Eqs. (15) in Eq. (14) 
and Eq. (7d) and solving for each term of the perturba- 
tion series one can immediately see that the zero order 
term in k is simply the DC Stark case of Eq. (8) repro- 
ducing the stark frequency uj e = vT+ A 2 , k = 1 term 
leads to the generation of new frequencies due to the 
possible two-particle interactions, i.e. self interaction re- 
sulting in second harmonic generation at 2tu e and 2uj as 
well as interactions due to coupling of electron and lattice 
at the difference wdfg = w e — w and sum wsfg = w e + oj 
frequencies. For k = 2, in addition to the frequency com- 
ponents already present for k = and 1 three-particle 
interactions give rise to frequencies at Sojij, \2uii ± ujj\ 
where i and j represent electron and lattice interchange- 
ably. This behavior is verified by Fourier transforming 
numerical solutions of Eqs. (7) and is shown in Fig. 4 
for the case when electron is driven close to the lattice 
resonance. In addition to the frequencies of two- and 
three-particle contributions analyzed above, higher order 



contributions of 4, 5 etc. particle interactions is evident. 
As mentioned in the previous section, there is a strong 
difference frequency component wdfg that is responsi- 
ble for the transport enhancement in the weak-coupling 
approximation (x *C 1). Finally, we want to point out 
that, as the resonance condition is approached, the ap- 
proximation e = fails. The qualitative effect of the field 
on the transport remains the same, but the fundamental 
and higher order electronic frequencies will depart from 
the values stated above. 
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FIG. 4: Frequency spectrum of the non degenerate non- 
adiabatic dimer normalized to the maximum close to res- 
onance. The different parameter values are \ — 0.005, 
A = 0.98. The initial conditions are the same as in Fig. 1. 

So far we have shown that the application of the elec- 
tric field results in transport enhancement through the 
generation of low- frequency wdfg component, where in 
via a two-particle process the electron is able to reso- 
nantly tunnel between the dimer sites, assisted by the 
interaction with a phonon. To underscore the role of the 
applied field in the generation of the transport enhance- 
ment let us consider the case of a degenerate dimer near 
resonance, i.e. when u e ~ uj = 1 (Fig. 5) without the ap- 
plied electric field. Since the delocalized initial condition 
corresponds to a stationary (non-oscillating) solution (as 
discussed in the previous section), we now consider an 
initially localized electron, i.e. po — 1, still in the weak- 
coupling regime (\ — 0.005). In this case (see Fig. 5) 
electron oscillations are supported at new frequencies re- 
sulting from the nonlinear mixing of the odd pairing of 
the two fundamental frequencies uj e and uj. These new 
frequencies are namely 2w — u) e , 2uj e — oj, 3oj e — 2uj along 
with the corresponding sum frequencies (not shown in 
Fig. 5). The frequencies corresponding to even pairing 
(e. g. wdfg) are completely absent. This is also evident 
from Eqs. (13), which, for the case of A = 0: 



(l + x y ) p = -xv r o 



x 2 y i i 
-Y\y 



u 



(16) 
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As before, the exact solution to Eq. (16) and Eq. (7d) 
is not possible. Following the perturbative analysis, the 
zeroth order term in p represents the equation of a sim- 
ple harmonic oscillator with no coupling to the lattice 
which has non zero solutions for not completely delocal- 
ized initial conditions (7*0 ^ 0). The first order correction 
represents the equation of two linearly coupled harmonic 
oscillators while the nonlinear mixing of frequencies can 
only be observed in the second and higher order terms. 
The newly generated frequencies result from the odd mix- 
ing of the two fundamental frequencies due to the cubic 
variation with respect to the oscillator variables. This 
situation should be contrasted with the case of non-zero 
electric field, presented by Eq. (14). There, the depen- 
dence on the lattice variables is quadratic, which means 
that while Eq. (16) and Eq. (7d) are symmetric under 
the transformations p — > —p and y — > — y Eq. (14) and 
Eq. (7d) are not. The fact that the potential remains 
centrosymmetric without the applied electric field pro- 
hibits the generation of the even harmonics. Application 
of the electric field results in the introduction of non- 
centrosymmetric component of the potential, hence al- 
lowing for even harmonics to be generated. 



Q. 
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FIG. 5: Frequency spectrum of the degenerate non-adiabatic 
dimer at resonance from Eqs. (7) for localized initial condition 
(po — 1). Other initial conditions are qo = ro = yo = yo = 0. 
The coupling parameter \ = 0.005. The two fundamental 
frequencies are denoted by uj e and oj. 

The impact of different initial conditions on electronic 
transport can also be understood from Eqs. (14) and(or) 
from Eqs. (lib), po = 1 initial condition corresponds 
to solutions in which the electron and lattice oscillations 
are completely out of phase, whereas ro = 1 corresponds 
to the symmetric mode or in phase oscillation of the two 
oscillating systems. Completely delocalized initial condi- 
tion as seen in the Stark case implies enhanced tunneling. 
These two effects coupled with the nonlinear interaction 
makes this initial condition more favorable for enhance- 
ment of transport. 



CONCLUSIONS 

In this paper we have presented a study of the response 
of the nonlinear quantum dimer to a constant externally 
applied electric field. Our studies differ from the previous 
ones (see for e.g. [21, 22]) in that we allow the lattice to 
evolve dynamically at its own time scale. Inclusion of the 
lattice degree of freedom results in a resonant enhance- 
ment in the transport of the electron. We performed 
numerical and perturbation series analysis in the limit of 
weak coupling of the electron with the lattice (x < 1) 
to provide a deeper insight into the ongoing processes. 
We find that application of the DC electric field induces 
a second-order nonlinear mixing of the electron (u) e ) and 
lattice (lu) frequencies facilitating electronic transport en- 
hancement near resonance. This mixing is manifested by 
a two-particle interaction resulting in a strong difference 
frequency component ujdfg — We~u, which corresponds 
to the oscillation of the slowly varying envelope of p (see 
Fig. 2). By symmetry arguments, lodfg arises due to the 
applied field and would not have been present, otherwise. 
The observed enhancement of the transport is similar to 
the phonon-assisted electron hopping, well-known for ex- 
tended systems [7] . Despite the classical treatment of the 
lattice, the subtraction of the lattice frequency u in the 
DFG process is akin to single phonon annihilation. This 
means that the electron absorbs one quantum of the lat- 
tice energy to complete the transfer to the next lattice 
site. 

We would like to further strengthen this last point by 
demonstrating the behavior of the system under the ap- 
plication of DC field with increased strength of coupling 
(Fig. 6). We show the variation in degree of transport 
with changing field (left), on a logarithmic scale and the 
contour plot of the different frequency components gen- 
erated. Apart from the already discussed resonant en- 
hancement w e ss lj, a sharp increase in 7^ is also seen at 
a higher field value (A rj 2.83), accompanied by a low- 
frequency contribution as depicted on the contour plot. 
This enhancement can be assigned to the field-induced 
three-particle process in which electron transfer to the 
next dimer site is assisted by the annihilation of the two 
quanta of the lattice vibration. Due to the higher or- 
der nonlinear process, the linewidth of the enhancement 
7d (A) is narrower than in the single phonon interaction. 
For stronger \ values, the enhanced transport is seen at a 
multitude of increasing field values, corresponding to the 
generation of low frequency components. In such cases, 
system quickly leaves perturbative regime and becomes 
chaotic, making it difficult to obtain a clear insight into 
the phenomenon. 

Our study should be relevant to systems which are lim- 
ited to a few lattice points, superlattices etc. but with 
the constraint that the system is initially not in a pure 
state. Even for extended systems such as the one treated 



7 



in the recently reported polaron study [14], a similar non 
linear mixing of frequencies might be observable. 
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FIG. 6: Spectral response for the probability difference as a 
function of the electric field (right). The value of the coupling 
X and lattice vibrational frequency uj was taken to be 0.05 and 
1.5 respectively. On the left is shown the degree of transport 
7d corresponding to that value of the field on a logarithmic 
scale. Other parameters are the same as in Fig. 1. 
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